






























































20-(continuation) 


v to quantum transport within the Liouville formulation. The second part con¬ 
ns with fluctuations of carrier velocities and energies both in stationary and 
transient regime, described by means of the correlation-function method. An 
analysis of the results, obtained through a Monte Carlo procedure, for cova¬ 
lent and polar materials yields a deep physical picture of the effect of the 
scattering mechanisms (phonon and carrier-carrier interactions) on the trans- 
| port properties. • >. „ y . f - f,! 
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1. INTRODUCTION 

The present report contains technical matter related to 
the research performed at the Department of Physics of the 
University of Modena for the Contract number 
DAJD45/83/C/0039 "Monte Carlo analysis of quantum transport 
and fluctuations in semiconductors". The subject is treated 
at three different levels. 

a) Part of the research has reached definite conclusions; 
in this case the material is well-organized and new results 
are presented. This the case of the work performed on 
fluctuations and on the effect of electron-electron (e-e) 
interaction on energy relaxation (Chapter 4). 

b) Part of the research is still under development and the 
related material is still in the form of models and 
proposals. This is the case of the research devoted to the 
attempt to extend the Monte Carlo method to the solution of 
quantum transport problems (Chapter 3). 

c) Part of the material contains a review of known results 
already present in the literature. We felt useful to perform 
this investigation and to present it in a unified way in 


orde 


in which our 









s report is 


formal quantum transport theory is reviewed in Chapter 2 
Different approaches to quantum transport such as the Maste 
Equation, the Generalized Langevin Equation, the Gree 
function method and the Wigner function method are describe 
and compared in their merits and shortcomings. 

Chapter 3 contains the description of the attempts mad 
to generalize Monte Carlo methods to quantum transpor 
within the Liouville formulation. A method is proposed t 
solve integro-differential equations of the type obtained i 
transport theories. The application of the method to th 
actual equations governing the motion of electrons in th 
presence of external and phonon fields is still unde 
development (Section 3.4). 

The research performed on fluctuations of carrie 
velocities and energies is described in Chapter 4. I 
particular Sect. 4.1 deals with the autocorrelation o 
velocity fluctuations, noise and diffusion in steady-stat 
conditions. while Sect. 4.2 describes the same Quantities i 












2. FORMAL THEORY OF QUANTUM TRANSPORT 


2.1 Introductory Remarks 

In the domain of ultrasmall structures (submicrometer) 
and ultrafast times (subpicosecond) simple arguments based 
on the uncertainty principle show that hot-electron 
transport in semiconductors needs a more exact approach than 
that offered by the semiclassical Boltzmann equation. Table 
2.1 reports the various levels of description of 
nonequilibrium statistical mechanics. Starting from the 
microscopic level, identified by the Liouville von Neumann 
equation for the density matrix of the whole system 
(electrons plus scattering centers), the kinetic level 
describing the time evolution of the electron system can be 
achieved. Then, by introducing several approximations, the 
semiclassical Boltzmann equation can be, for example, 
obtained. 

Theoretical efforts towards more exact quantum approaches 
have indeed predicted new phenomena /2.1/. Among these, the 
intra-collisional field effect (ICFE) and the collisional 











The former accounts for the fact that the scatten 
occurs between electron states in the presence of 
electric field. The latter is correlated with the fini 
lifetime of the electron states as a result of collisions 
The main consequence of both effects is in the expectati 
of a smoothing in the sharp peaks of the scattering rat 
occurring at threshold energies, as depicted in Table 2.2 
An illustrative description of the above effects, 
compared with the semiclassical Boltzmann picture, 
offered in Figs. (2.1-3). 

Within the semiclassical Boltzmann picture (see Fi 
(2.1)), an electron performs a classical trajectory in re 
space between successive collisions, which are treated 
point-like events in space and time. Transitions betwe 
initial k and final k' wavevectors are also point like, a 
the well defined correspondence between k and £(k) enabl 
a straightforward calculation of the energy involved. As 
result of this picture, the transition rate P(k,k') 
proportional to an energy conserving delta function. 

In presence of ICFE the concept of "duration of 
collision" 72.2/ can be usefully introduced (see Fig. 2.2) 






- — c (q being the transferred wavevector |k-k‘ | and 

the duration of the collision) can be defined /2.3/. In 
this way, the carrier can gain or loose energy from the 
field during the collision. Thus, owing to the 

interference between the field and the scattering mechanisms 
and the possible gain of energy during the collision, a 
scattering event can occurr which in the point-like 
collision model was forbidden by the energy conservation. 
As a result of this picture, the original energy conserving 
delta function is transformed in a more complicate 
expression (see Fig. 2.2), which can be analyzed in terms of 
Fresnel integrals /2.3/. 

Collisional broadening is linked to the manybody nature 
of the phenomena. The quasi-particle approximation must be 
released, and a complex self-energy 2L (^, e ) has to be 
introduced (see Fig. 2.3). As well known /2.4/, the real 
part of renormalizes the energy of the electron state, 

while the immaginary part, through the optical theorem, 
gives the scattering rate. Energy and wavevector are now 
two independent variables and a spectral density, S(k,€ ), 
giving the probability of finding the electron with (k,<£ ) 


must be introduced. 


As a result of this picture, the 


















original energy conserving delta function is substituted by 
Lorentzian structures. 

Different lines of approach can be pursued to establish 
quantum kinetic equations able to describe the above 
phenomena /2.5/. The aim of this chapter is to present the 
main results which, to the authors opinion, have recently 
appeared in the literature in an attempt to evidence merits 
and shortcomings of different formal theories. Our hope is 
to obtain indications on approaches which are most 
appropriate for numerical solutions, having in mind the 
possibility of extending to a quantum mechanical framework 
the traditional Monte Carlo technique /2.6/. 

The chapter is organized as follows. Section 2.2 will 
specify the physical model. Sections 2.3 to 2.6 will survey 
the four main lines of development of the formal theory. 
Applications and results will be reported in Section 2.7. 




2.2 The Model 
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We concentrate on a simple case, yet without any loss of 
generality, which corresponds to an ensemble of independent 
electrons. These electrons are under the influence of an 
applied electric field E which, if not otherwise stated, is 
taken as homogeneous in space and independent of time, and 
they interact with a phonon bath. Thus the total 
Hamiltonian of the system is partioned as; 

H = H* + H f + H P * H t .p (2.i) 

where H g refers to electrons Hj_ to the electric field, H p 
to the phonons and H e ^ to the electron-phonon interaction. 

As the solution of a transport problem relyes on the 
calculation of average values of the observables of interest 
(e.g. density current, concentration, energy, etc.) a 
prescription for such a calculation is needed. To this end 
different approaches can be used and, in the following, the 
leading ones will be surveyed. 
















2.3 The Generalized Master Equation 


2.3.1 The Generalized Master Equation 

Within a Schroedinger picture, the density matrix (t) 
of the physical system of interest is introduced. The 
average value at time t of an observable 0, which does not 
depend explicitely on time, is given by: 

< 0> fc = T* (2.2) 

The density matrix obeys the Liouville von Neumann equation: 

^ (t) = “ —■ -L ? (2.3) 

L = £h , (square brackets denote commutator) being the 
Liouvillian superoperator associated with H. 

Following Nakaijma and Zwanzig /2.7/, irrelevant 
information about the electron subsystem can be eliminated 
through time independent projection operators. In this way, 
in place of Eq.(2.3), one can write a kinetic equation for 


density 


(gel 









uation), which can be analyzed perturbatively in terms of 
e applied field and of the electron phonon interaction, 
ing this procedure Pottier and Calecki /2.8-10/ have 
troduced projection operators that allow a factorization 
the density matrix for the coupled electron phonon system 
d they have obtained an interesting quantum kinetic 
uation. To restore the traslational invariance of the 
ectron Hamiltonian in presence of a uniform electric 
eld, they used a vector potential gauge; 

+ H f = JL-[ p - (2 * 4a) 



A (t) 


-J fc E(b') ib' 


(2.4b) 



ere m is the electron mass, p the momentum operator and 
t) the vector potential operator. Furthermore, 
atistical electron operators are introduced which at t=0 
e given by; 


ynskjo) = 



IRXfcl 


(2.5) 










'here time t=0 corresponds to the switching on of the field, 
ind thejK ^ states are the accelerated plane waves. If at t 
' Y is evoluted by the total Hamiltonian then, to the lowest 
irder (second) in the electron-phonon interaction, the 
quation of motion for Y(K,K;t) is obtained: 

. _ _ t _ _ 

y(R,K;t) = I At'Z [P(K,Z;t,t')Y(K,K;k‘) - 

>o k' 


- * F(E;tJ ,2.6.) 


here the probabilities P(K,K';t,t') entering the integral 
ollision term are given by: 


P(F,5';t,n = f, (?* I I I** 

» \C ACn) - 

- € (K- 4 AOT) -/y*<«JJ} (2., 


1 


1 (-1) corresponding to phonon emissii 


(absorption) 







processes. F(K;t) represents a stochastic (Langevin) force 
which is given by: 

FfM) = Z K . 

[CK'IY,/S>Ya',i<,o) tiplif'Jt'fed'-gAtt'ij-efi-t AW)]- 
.<1 |^|R'>Y \*cU'(e(t;'-l A(PI)-e(R-< AU'))]]J 

(2.6c) 


In contrast with the collision term, F(K;t) does not 
involve the projections on the JK states, but only 
off-diagonal electron operators. In Eqs.(2.6) C(q) 
characterizes the nature and strength of the electron phonon 
coupling, ^j= expticj • ’r), is the angular frequency of 


the phonon 
population, 


£(k) the eigenvalues of the zero field 


Hamiltonian, and'tiK the canonical momentum given by: 


fiK = -hk + eA( t) 


( 2.6d) 


The electron distribution function f(K,t) in I K ) states is 
















conveniently expressed in the Heisenberg picture and results 


to be: 

« T* (2.7) 

where the trace operator acts in the whole electron phonon 
space. 

By assuming that the initial density matrix (0) is 
factorized in the phonon and electron parts as: 

? (0) = ?ph-fel CO > |2 - 7a) 

one gets: 

Tx. [f(o) r(K;t)J =0 <2-7b) 

So the Langevin force does not contribute to the evolution 
equation of the electron distribution function. Therefore, 
f (~K; t) obeys the high field retarded transport equation: 






(2.7c) 


As merits of Eqs.(2.6): (i) a Boltzmann Langevin equat 
able to describe average quantities as well as fluctuate 
around average is obtained from first principles. (ii) 
explicit microscopic expression of the random Langevin foi 
is provided. (iii) By including the external field in 
basis state, the intracollisional field effect 
automatically accounted for. (iv) The choice of the vec 
potential gauge enables to avoid the use of the A 













2.3.2 Comparison Between The Generalized Master Equation And 


The Boltzmann Equation 

To gain some physical insight, we shall compare Eq.(2.7c) 
with the familiar Boltzmann equation: 

oictt.t) ,i DjcK.t) = Z., fpotlviHKV- 
~Ti~ * l u L 

-r(tX) (2 .7d) 


To recover Eq.(2.7d) from Eq.(2.7c) the canonical K-vector 
has to be substituted by its expression as a function of the 
crystal wavevector (see Eq.(2.6d)). In doing so it is 
easily verified: 


1 J&t) (2.7., 

ot ot k oe 

I t /** 

Then, by taking the ansatz of completed collisions , J I 
, in the r.h.s. of Eq.(2.7c), and neglecting the effect of 
the field during the collision, the collision term 
P(K,K';t,t') recovers the usual markovian form in the 


carrier wavevector representation P(k,k'). 







2.4 The Generalized Langevin Equation 


2.4.1 The Generalized Langevin Equation 

Within the Heisenberg picture, the average value o 
Eq.(2.2) writes; 


< 0> t = T* [f(0) 0 H (t )j (2.8) 

where the time evolution of the operator 0 ( (t) (here th 
subscript H indicates the Heisenberg picture) is given b 
the Heisenberg equation of motion: 


6 U = it o H (t> 


Following the original idea of Mori /2.11/, Zubarev /2.12/ 
and Grabert /2.13/ , a decomposition procedure which make 
use of time-dependent projection operators can b< 
introduced. In this way a generalized Langevin equation fo 
the thermodinamically relevant variables (macrovariables) o 
the problem can be written in place of Eq.(2.7), which cai 


be analyzed perturbatively. Following this procedure, Ferr; 












coworkers /2.11-13/ 


have 


considered the set 


:rovariables: 


[ ] ~ j j Mp , H^.p] (2.10) 

■re P is the total momentum of the electronic system, 
> number operator of electrons and the ot 
:rovariables are the same as in Eq.(2.1). Then, 
lation of motion for the macrovariables is found to be: 


P(t) = exp (i Lt) 'tf(L) F(o) + 


+ fds 1 <up (cLs'} i V(s') L[j-7T fs'jJ 6 (s,t) P(o) + 

J o 

r- * ft , 

U-Vtoil 6(0,t) P(0) - ds. t,f> (iLs'l r(s') 6 (s,t) 
,[d~tr(t)]F(o) 


e IT(t) is a time dependent projection operator which 
ined in terms of the chosen macrovariables and of their - 







average values, and the operator G(s,t) is a two time Greei 
function given by: 

Q u.t) = r_ t*t> [ <•' j‘ l I> wul] l2 - 121 

T being the Dyson time ordering operator from left ti 
right. 

The physical meaning of the terms on the r.h.s. o: 
Eq.(2.11) can be identified in the following way. The first 
term gives the collisionless motion. The second term can b< 
divided into two parts, the former characterizing th< 
collisions and the latter leading to a fluctuating forc< 
induced by the electron-phonon interaction and th< 
nonequilibrium nature of the system. The third and fourth 
terms are two additional fluctuating forces, the forme: 
being Mori like in nature and the latter characterizing th< 
fluctuations induced by the rate of change of th< 
macrovariables during the transient regime. The averagi 
motion of the macrovariables is determi >ed by the equation: 


<P(t)y = < ex-p(iLl) irffc) p(0))> 






ds'ixp (CLs!) Tf (s') i L [j-irCs'jjGW) P(o))(2. 13 a) 


As merits of Eq.(2.12): (i) A generalized Langev 
equation valid at high electric fields is obtained fr 
first principles. (ii) Separate descriptions of the moti 
in absence of collisions, with collisions, and wi 
fluctuations are evidenced. (iii) Fluctuations und> 
transient conditions are accounted for. (iv) Generalize 
balance equations can be obtained. 

As shortcomings: (i) The choice of the pertinei 
macrovariables is critical but no fixed rules govern th 
choice, rather one must be led in doing that by physic. 





2.4.2 Comparison Between The Generalized Langevin Equation 
And The Classical Langevin Equation 

To get some physical insight, a comparison of Eq.(2.11) 
with the classical Langevin equation and with the equation 
derived from first principles by Mori /2.11/, under linear 
response conditions, is here reported. 

The Langevin equation, phenomenologically derived in 

presence of an electric field, takes the form; 

p (b) = ^ fCt) +K(t) (2.13b) 

where "p is the classical momentum of the particle, a 
relaxation time accounting for friction effects, and R(t) 
the stochastic force associated with the rapid fluctuating 
part of the carrier motion. A first principle justification 
of Eq.(2.13b) under linear response conditions in the 
external field was given by Mori /2.11/. By introducing a 
time independent projection operator IT Mori's equation for 
the quantum momentum P writes: 


P U) = 4rp(iLt)'i7 P(o) 








( As arpQUt-s) JTT^L trpW-nJLsJU-TfJPCo) 

Jo 


+ £xp[('{-7f Ji t] C-i-Tf) P(o) (2.13c) 

where the three terms on the r.h.s. of Eq.(2.13c) are the 
quantum analogous of those in the r.h.s. of Eq.(2.13b). It 
is worth mentioning that in the Mori equation the friction 
term is described by a quite general retarded kernel. 
Furthermore, first-principle justification of the stochastic 
force, postulated by Langevin, is provided for the first 
time. 

Under far from equilibrium conditions, the first three 
terms on the r.h.s. of Eq.(2.11) are the analogous of those 
in the r.h.s. of Eq.(2.l3.c). The nonlinear condition is 
however responsible for the need of a time dependent 
projection operator and for the fourth term in Eq.(2.11). 
This term is related to the transient which goes from the 
switching on of the field, assumed to be step-like, to the 
stationary far-from-equilibrium condition. 








2.5 The Green Function Method 


2.5.1 The Green Function Method 
Within a second quantized Heisenberg picture, with creat 
; (r,t) and annihilation l | / (r,t) field operators, 
average value of Eq.(2.2) writes: 


<o> t = 


Tn. u’CtJOtlCU] ^ 

T * (Ho-/•*'<)!] 


where the average is performed on the grand canoni 
equilibrium ensemble (i.e. H 0 = H - ^ = 1/(K (S T ), 

being the Boltzmann constant and T. the bath temperature, 
is the chemical potential). All the dependence of 
external field is esplicitated in the evolution operator: 

\ct) = % y 1 #,*') f(*U') h f ]\ {2 - 15) 

- a? 

T being the time order operator from right to lef 

Following Kadanoff and Baym /2.17/, a set of real time Gr 


functions is introduced as: 











J = - 1/ < l YClfb) (2.16a) 

Gf = i <f T (l,t) (2.16b) 

where > < indicate t ><t', respectively. The average in 
Eq.(2.14) can be calculated in terms of and G. In 
particular for the density current it is: 



C^-%) ]. 


tl-n,' 


On this basis, Wilkins and coworkers /2.18-20/ have 
developed a Generalized Kadanoff-Baym (GKB) formalism which 
allows to treat the external field exactly. Conceptually 
the calculations split into two steps. First, the motion in 
absence of collisions, but in presence of an electric field 
is treated exactly. This leads to the introduction of field 
dependent Green functions G^ tv : 









+ H c C, 


(2.18a) 







where ^ refers to retarded and advanced Green functions 
for the free electrons. Second, the description of the 
distribution of the particles in presence of collisions is 
given in terms of the field dependent Green functions. These 
functions satisfy the GKB basic equations: 




. Wp 




[(«V, C?] = * [ 2 *. 41 - 

_ d ( 2.160 

To have a closed set of equations one needs to know the 
electron self energy as a functional of the Green 

function G. To complete the notation let us remind that the 
following identities are valid for G (and ): 


sn,n = j 4 < 5 * ( 


(2.19a) 








•) = 0 (t'-t) [<z*(t,t') - 


(2.19b) 


Gf n (k,k') = d(t-t')[4 > (b / k , )-<9 < (b / b')~\ (2.19c) 

Q(t-t') is the unit step function and j, j indicates the 
anticommutator. Furthermore, as all G's are double time, 
double space, and real time Green functions, Eqs.(2.18) 
should be interpreted as matrix equations: integrals over 
intermediate position and time variables are implied 
throughout. In Eqs.(2.18c) all external forces are denoted 
by the operator Hp_ , whereas all scattering processes 
(impurities, phonons, etc..) are collected in the self 
energy. Eq.(2.19c) is a generalized quantum kinetic 
equation. The physical meaning of different terms is the 
following. The left hand side gives the collisionless 
evolution of the Green function. The right hand side 
represents a generalized collision integral: the commutators 
describe the kinetic effects of the interaction 
(renormalization of the electron eigenvalues due to the 
presence of field and interaction); the anticommutators 


be the dynamical ef 









obt, 





2.6 The Wigner Representati 


2.6.1 The Wigner Representation 

The one-particle Wigner function f^(p,r;t) /2.22-24/ is 
defined as the Weyl transform of the density matrix (t ] 
divided by h^ . In the case of a system in a pure state, f^ 
has the well known form given by Wigner: 

j(p,*;t) = ~ jdv 

P *) ffi '} 12 ' 201 

where l f , (r,t) is the wavefunction in configuration space. 

The average value of Eq.(2.2) becomes: 

S0> t = JJ a (p,*) jp e/c (2.21) 

where the function o(p,r) is given by the Weyl 
correspondence rule: 

& (p,*')- Tn I o fj da, ± ( (%-$) -u v- 



(p-fi)v-)]] 


( 2 . 22 ) 






















m (R 


G,(P,Z;k,T) = 


R-S T-£ ) (2.24b) 

£ ' * ' 2 2 J 


it 


is: 



f(PR,TJ = ^ jj„ tU*.?;".?) 


Thus the kinetic equation (2.18) is appropriate for 

introducing all relevant quantum effects into f(R,P,T). 
Following this procedure Barker and coworkers /2.26-28/ have 
analyzed collisionless (ballistic) transport for particles 
interacting with model potentials. 

As merits of the Wigner representation: (i) One can take 
advantage of the correspondence rule to introduce concepts, 
like Wigner trajectories /2.29/, which visualize quantum 
effects in a classical framework. 

As shortcomings: (i) The interpretation of the Wigner 
distribution should be made carefully as, by not being 
necessarily positive definite, it has not a probability 
meaning. In this respect a generalized Wigner distribution 
f(l*,P; oj , T) = -i (R,P; CO,T), which allows for an 

autocorrelation function interpretation, should be more 








onveniently introduced /2.17,30,31/ 
dentity between the generalized W 
Ireen function formalisms is clearly e 
ii) The stationary Wigner functioi 
rom the transport equation, un 
onditions are restricted /2.24,32/. 
n eigenvalue equation reviewed in Ref 







2.6.2 The Wigner Function For The One-Dimensional Potential 












Y = % (Vo-*) 


(2.25c) 


and A is an arbitrary normalization constant. 

By assuming that the system is prepared in the state 
represented by Eq.(2.25b) at time t=0, the corresponding 
Wigner distribution can be obtained from Eq.(2.20) applied 
to the one dimensional case. The result is; 




Vm aea 


_ 4 [fll Z 


<Lxp 




4 r k 2 


[l‘f ‘‘I'ripif 






VI# 1 


ij <£j>yO 


( 2.2 5d) 


and 


= 4 JA ] ‘ j cos (if?- ■?*«,) _ r ^ 

‘ fc [x‘f(ip-k)*]lr 


cos (iP? + i’H) 


r 


- - Sc'/H f 2 P° l _ f \ x 

lx'*(-i-P>l<F]Lr ' ft 1j 








k (j UfjJl 1L l _ _ si/h I *P? t. x 

[r'*( e £-VI# (%->') k 

_ (Zi£zi2 i 

[r^Cf^n %’(*>!') 1 

( 2.2 5e ) 

Since the particle is in an energy eigenstate, Wigner 
ajectories are given by equi-Wigner curves. Fig. 2.4 
ows the results. The classical trajectories for the 
esent case (see Fig.2.4(b)) consist simply of straight 
rizontal lines. An interesting feature exhibited by 
gner trajectories shown here (see Fig.2.4(c)) is that a 
n negligible portion of the trajectories penetrate through 
e potential step. In general, the larger the energy of 
e space-phase point approaching the potential step, the 
eper i.s the penetration . Very roughly, the average 
netration can be estimated by looking at the trajectories 


<:f y<° 



d their relative weights given by f. This approach thus 









om Fig.2.4(c) some trajectores are associated with 
gative weights. This is due to the fact that the Wigner 
stribution function is not positive definite. The Wigner 
stribution function is essentially an auxiliary 
thematical function introduced for convenience of 
scussion. For the sake of completeness. Figs.2.5-7 report 
e full Wigner function in the classical phase space for 
s different particle energies there indicated. 



fir' i 
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2.7 Applications 

Real cases can be classified within two limiting regimes 
of motion: the time reversible collisionless case and the 
irreversible stationary case when a given particle undergoes 
many completed collisions. The formal theories outlined in 
the previous section cannot be applied "sic et simpliciter" 
to the latter case. Indeed to identify a given number of 
completed collisions a coarse-graining of time is required. 
This problem concerns the reduction of the original 
non-Markovian generalized collision integral to a Markovian 
form. This procedure, which is at the basis of 
irreversibility, is however associated with a loss of 
information which is hard to evaluate in terms of physical 
approximations. This is the actual bottleneck of quantum 
transport theory. 

Using the Green function method as a guide line, an 
interesting attempt to account for quantum effects within a 
standard Monte Carlo procedure has been proposed by Hess and 
coworkers /2.34-36/. A self-energy model for an electron 
nonpolar optical-phonon scattering mechanism is introduced 








Zte)= f}dc 


(2.26) 


£ — £ —'hu>ap—2L(€-ku , Op)+ l 'S 

where g is the coupling constant given by 'TidV ( 2 t-0^) (h 
being the Planck constant divided by 2ff , D is the optical 
phonon deformation constant, ^ the density of the material) 
and |3 (£. ) the electron density of states. 

Eq.(2.26) is solved numerically using a density of states 
obtained from pseudopotential calculations. The optic 
theorem allows to define the scattering rate for t 
electron phonon interaction as: 



= - rf ) I/w Z. (e) = - (%) T 1 (2.27) 


which should, in principle, improve the Golden rule result 
Indeed, self-energy corrections are found to reduce t 
scattering rate. Furthermore, the delta functions appeari 
in the energy conservation terms of the Golden rule are 
replaced by Lorentzian structures of the type; 


S (<£<,-£<. ) -► 


i * n) 

ir u L - 6^ +(r c +r 5 f 


(2.28) 


m - 
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states, and €. ■ / are the quasi-particle energies shifted 
by the real part of the self energy. Within a Monte Carlo 
procedure, Eq.(2.27) can be used for the determination of 
the free flight (self-energy correction) and Eq.(2.28) for 
the determination of initial and final states (collisional 
broadening). The latter step is a delicate one since it is 
not clear how it could be included consistently into 
numerical calculations. This point needs a more careful 
analysis than that presently available. 

Several authors have recently speculated about the field 
dependence of the transition rates in high field quantum 
transport. Within a Stark ladder representation Sawaki 
/2.37/ has analyzed the effect of an electric field in terms 
of the self energy due to higher order terms in the elecron 
phonon interaction. Herbert and coworkers /2.38,39/ have 
used a time dependent perturbation theory modeling the 
scattering between Airy functions. Analogous calculations 
have been performed by Ziep and Keiper /2.40/. Marsh and 
Inkson /2.41/ have attempted to generalize the above results 
by employing WKB wavefunctions in place of Airy functions. 
Overall, no definite conclusions have been reached so far, 


and different approaches are still a matter of discussion 









2.8 Conclusions 


Recent developments on quantum transport in the 
hot-electron field has been surveyed. By comparing 
different formal approaches, the Green function method 
emerges as a more powerful technique since it can take 
advantage of a diagrammatic perturbation expansion for the 
interaction of carriers with the scattering centers. The 
comparison between theory and experiments /2.34, 35, 
43-45/is still far from being satisfactory, in view of the 
difficulties to estimate the reliability of the 
approximations introduced and to find crucial experiments on 


the subject. 
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OWING TO SCATTERING PROCESSES THE ELECTRON STATE HAS A 
FINITE LIFETIME. THEREFORE THE ENERGY WAVEVECTOR RELATIONSHIP 
IS DESCRIBED BY A SPECTRAL DENSITY S(1T, £ ). S(iT, £. ) IS THE 
PROBABILITY OF FINDING THE ELECTRON WITH GIVEN ("k, 6 ) . 



SHARP PEAKS OCCURRING AT THRESHOLD ENERGIES ARE SMOOTHED 




TAB. 2.2 











FIGURE CAPTIONS 


Fig. 2.1: Schematic representation of the semicl 
Boltzmann picture. 

Fig. 2.2: Schematic representation of the Intra-Co 
Field Effect. 


Fig. 2.3: Schematic representation of Collision Broade 


Fig. 2.4: Representation of the Wigner trajectories 
one-dimension potential step model (arbitrary uni 
used) . 

(a) Profile of the potential. 

(b) Classical trajectory for the case of a particle energy 0 
lower than the potential energy of the barrier V u . 

(c) Wigner trajectories for the case 0 = 0.5, V 0 = 1. 

Fig. 2.5: Plot of the Wigner function for the one dimension 
potential step model in phase-space. £=0.1; V 0 = 1 in 










































FTG. 2.4 




































3. 


MONTE CARLO SOLUTIONS OF TRANSPORT-LIKE EQUATIONS 


3.1 Assumptions and limits of the traditional Monte Carlo 
procedure 

The Monte Carlo method, as applied to charge transport in 
semiconductors, consists of a simulation of the motion of 
one or more electrons inside a crystal, subject to the 
action of external forces due to applied electric and 
magnetic fields and of given scattering mechanisms. The 
duration of the carrier free flight (i.e., the time between 
two successive collisions) and the scattering events 
involved in the simulation are selected stochastically in 
accordance with some given probabilities describing the 
microscopic processes. 

As a consequence, any Monte Carlo method relies on the 
generation of a sequence of random numbers with given 
probability distribution. 

When the purpose of the analysis is the investigation of a 
steady-state homogeneous phenomenon, it is sufficient in 
general to simulate the motion of one single electron: from 
ergodicity we may assume that a sufficiently long path of 


if OJ 


























the physical systems of interest have been analysed and 
understood within the limits of validity of the 
semiclassical theory of transport. 

Therefore, a need is now felt to overcome these limits and 
try to extend the method of Monte Carlo simulation to 
include quantum effects not accounted for in the traditional 
approach. 

Such a need is further justified, on practical grounds, by 
the present state of VLSI technology which now allows the 
fabrication of devices so small that quantum effects may 
become appreciable. 

Quantum effects are a consequence of small dimensions not 
only because of the well known size effects, but also 
because the applied voltages (which cannot be decreased 
below certain limits) at such distances imply very high 
fields, at which quantum phenomena may become important 

In the traditional semiclassical approach, collisions are 
treated as pointlike events in both space and time. In 
fact, a fully quantum mechanical treatment should account 
for the interference between applied and scattering fields 
(intracollisional field effects), the finite lifetime of the 


electron 


states 


(collisional broadening), 


and 


the 
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possibility of multiple collisions. 

Within this picture, energy conservation is not a feature 
of each of each scattering event, and the simple one to one 
correspondence between energy and wavevector will be 
substituted by a probability distribution ( £. ). 

Furthermore, at any new collision the assumption is made of 
perfect energy conservation of the previous collision, which 
is not exact when two scattering events happen quite near to 
each other. Collision broadening is thus neglected. 

Another feature of the traditional approach which we do not 
feel happy about is the use of first order perturbation 
theory, that is the use of the Fermi golden rule, for the 
transition probabilities. 

The situation then calls for a search of new techniques 
which can allow us to treat transport problems with a fully 


quantum approach, free from the above limitations. 


















the series in Eq.(3.2) converges/3.2/. 


It may be interesting to put Eq.(3.2) in the form 


(I-A)x = f (3.4) 


or 


x 


(I-A)"' f 


(3.5) 


In this way we recognize the series expansion in Eq.(3.2) a 
the harmonic series expansion of the resolvant operator i 
Eq.(3.5) : 


(I-A) A = 1+A+A !L +... (3.6) 

A Monte Carlo procedure has been developed by Von Neuman 
and Ulam to solve the system in Eq.(3.1), which proceeds a 
follows/3.2/. 

Let us consider a set of n "states" or "boxes" labelled 1 
2, ...n, and let a "simulative particle" be positioned a 
the beginning of the simulation in box i. Let then simulat 
a terminatina random walk with fixed but arbitrar 










transition probabilities 


p *j v i* j = lf • • - n 


from state i to state j. p.. must satisfy the following 


p ^.0 (>0 if A 4 j/0) 


2-4 Pa< 1 i=l, ... n 


For each state i the random walk has a terminating 
probability given by 




ftny random walk ^ may be represented by the set of the 
successive visited states 


"P = ( i. , i.i . , i. ) 




given random walk represented by 
>ute the quantity 


Eq.(3.10) let 


ipression 
walk has 
pectation 
is given 


P< V/io) 


X(^ ) = V( Y )• 


Ht 


(3.11) 


V(^\ ) 


Auil A if 


b. . P . jo 

1 to tt / / .’ ' 


V( -p = l if ) (3.13) 

= (i 0 ) in the last equation means that the 
terminated before any transition took place, 
value of X( ) for random walks starting 
by 

X(f/i„ ) =£, P«7i„ )X(C/i„ ) (3.14) 

< 7(0 

is Lhe probability of having the random walk f 







< Xte/*0> - Z_4 £-1 t-A Pfo/'o) 

k=o t-t. tfc. 
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H 

Z_^ 2L4 •" 4 _i 

t=o it tk_ 


f>. . b . • •• b . 

1 wii iaa I„ J 


A*>l 
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f + ^ ^ioiL fa + < r-r i A«n f u -f 

ti J tiC£ * 

£ + Af+ A f ‘ - ^io 


ere we have taken into account Eq.(3.2). 
erefore the following Monte Carlo procedure 
order to solve the system in Eg.(3.1): a gi 


chosen and many random walks are performed 








the state i, 


, according to the ti 


probabilities given in Eqs.(3.7-9) 
the quantity X( ^ )in Eq.(3.11) is evaluated. Its average 



value over the random walks is a corr 


















3.3 Generalization _of_ the _ procedure _to_an 

integro-differential equation 

Let us consider an integro-differential equation of the type 


b 







where v is a constant and H(x,x',t) is a known function; 
f(x,t) is the unknown function defined in the interval 
a^xi'b. 

The interest in considering this equation is based on th< 
fact that transport equations have often this form. Ii 
particular the Boltzman Equation can be reduced to the fori 
of Eq.(3.16) with H closely related to the scatterin< 
probabilities and v proportional to the applied electrii 
field. 

In what follows we shall describe a numerical proceduri 
which generalizes the method shown in the previous sectioi 
to find a steady state solution, if it exists, of ar 
equation of the type of Eq.(3.16). 



Eq . (3.16 ) we may obtain 
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b 

f(*. htik)-f(x,h) - 1 rifc 2^ + Jl4(*.md:(fa)<h} (3 . 17) 

a 


or 


b 

^x,fc4t) = fu-i*, b)+ Ju(xXt)l!L px't)A' 

o) 


(3.18) 


where x=vdt. 

Let us consider a mesh of the interval (a,b) and let f^(t) 
be the value of f(x,t) in a representative value of x in the 
i-th interval. If x is in cell i, let x-^x be in cell i'. 
Then our equation reduces to 


f./UHA,. fj(t) 









(3.19) 







Eg.(3.16) has an 


asymptotic solution independent of t, 


large enough t Eg.(3.19) reduces to 


(3.21) 


3 eguation is not exactly of the type desired (Eg.(3.21)) 
:e f-i is not known. However an iterative procedure can 
attempted which starts with a trial function f ; If f 
inserted in the first term of the r.h.s. of Eg.(3.21), 
with the Monte Carlo procedure indicated above 
tion f*^ can be found, solution of the equation 


(l) Jo) 

f; = f;i + 


(3.22) 


process can be repeated with f as known term, to fine 
next iteration f ' , and so on until convergence n 






tested the above procedure for the case of 


-(Ax' 1 )/ L 

fT t (& fcxV-y-trW*’)' 

= Vlr (x'6 b -t \ ) 


here x is variable in the whole real axis. In this case 
he equation has the solution, nomalized to unity, 


f(x,t) 



(3.24) 


fith asymptotic solution given by 


f(x,t) 



-X 

e 


(3.25) 


rhe Monte Carlo procedure described above has been applied 
>nd the results are reported in Fig.3.1. Here the initial 
:rial functionand the function obtained after a large number 
>f random walks are shown together with the exact solution 


»f the equation. 








Let us consider an ensemble 


the presence of a homogeneous static electric field E ai 
perturbation potential given by a superposition of p 
waves with different wave vectors and equal frequency. 

w/ l\. y L * -ifr-u. t)l 

V^|t) “ Afj & + /A(j & j ( 3 . 26 ) 

This potential may be considered a semiclass: 
representation of an optical phonon field. 

In what follows we shall adopt a density matrix formal: 
The one-electron wave functions can be expanded in a sei 
of Bloch states as 

. c- _ ti) . t k r 

Y o.t) ' 2-1 C (l,i) M.(f) 6- < 3 - 271 

It 


where i indicates the i-th electron. 
The density matrix ^ is defined as 
























Tf Zj b (C,t)/M fe ff)e/ 

It 


k>'V £ .t) = C w ft fe)6 






The density matrix is now written as 


M , l)f 

(kt't) •’ Z<L b'fttlb /E,'t) - ^£,fe',fc)& <3.331 


By substituting ^ as given by this expression in Eq.(3.30) 


, / -i(uh-C0*)i\ 

aW)t 



we obtain 









(3.34) 


-t / U\.")b 

a- H(t“b'b) 


. f r itr ,t'r 

H(bUb)* Jjf/Ut(f)& H/w^e- = 

/ft. * . -itr jWtfc 'br -/<<£, t| 

* ( jiYMjxie' e- H/Wt'Crte- ^ 

-i (M.-hV)b -<(uU-k)J)k. 

• b- = ^i(iL'b) e- 


Now ^(j(k,k',t) are the matrix elements of II between the 
dependent Bloch states. For simplicity these will 
approximated, in what follows, by plane waves, 
neglecting the overlap integrals in the matrix elements. 
By inserting Eq.(3.35) into Eg.(3.34), we have 


L / '?aCkt!b)\ -I 

* k -rr— «- + W< 

















^ jgl(U'fc) 7 ‘'“ a 4"k'b) e ^ ‘”'".3.3 

- cl [ik h)t/ 7 tfcOfifye,' 


or 


SCLto) , , . . _ f 

^ -^r _+ ' k/Wt ' £o ‘ ,)a4Lt) - % KfaWty 
, \ .. * /.0 *=■ / 


- a tot) 


The total hamiltonian is given by the sum of three terms: 
*K< 0 = unperturbed hamiltonian for an electron in the perfei 
crystal; 

= hamiltonian due to the presence of the external fieli 
here assumed uniform and constant; 

= hamiltonian due to the interaction with the oscilla 

(phonon) field given by Eg.(3.26). 


For the matrix elements we have 
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In order to consider V as representative of the phonon field 








I 




(3.41) 


where ^ is the density of the crystal, V its volume, th 
deformation potential coupling constant, and N^ the phono 
occupation number. 

Furthermore, in the second term, corresponding to th 
emission of phonons must be substituted by + 1, i 
order to take into account spontaneous emission. Therefore 






~ {to? 


-i [U>t> -iO c -lJ 0 )L 


By substituting the matrix elements calculated above inti 
Eq.(3.37) for the dynamics of the density matrics, we obtain 


= -k A CLfLL'y) 







-+ ^ Ji M - /— i-'fo-UL-H) 
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]/<2 V / 






4<i & 
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here k = eE/ff and 

his equation can be transformed into 


t) « CL Sk, fe- < th, t-st) + 

-I* U- (Vow)tSfc. a Aft) 

....{£ '...W; 

ViktJAkf\J i ' 

-($]■■■■- a(u"t) 


ich can be reduced to the same form of Eq.(3.1). 
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FIGURE CAPTIONS 


Fig.3.1;Exact solution (dotted line) and numerical so 
(continuous line) as obtained following the p 
dure described in Sect.3.3, for uniform initi 
trial function (dashed line). 
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obtained 
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(Sect.4.2) phenomena, and also to analyze the different 
contributions to the diffusivity due to the different 
physical sources of fluctuations which arise in the presence 
of an applied electric field. 

Results have been obtained with a Monte Carlo procedure 
for covalent (silicon) and polar (gallium arsenide) 
materials. 

Furthermore, Monte Carlo calculations of hot-electron 
phenomena usually do not take into account interaction among 
carriers. On the other hand , a strong dependence of energy 
relaxation time on carrier concentration is observed at 
rather low carrier concentrations and it is thought to be due 
to carrier-carrier interaction. 

An analysis of the influence of e-e scattering on the 
energy relaxation time is presented (Sect.4.3.1) through the 
energy autocorrelation function obtained from a Monte Carlo 
procedure which includes e-e interaction 
(Sect.4.3.2). Results will be shown for p-Ge and compared 
with experimental data (Sect.4.3.3) . 

4.1 Autocorrelation of velocity fluctuations, noise , and 
diffusion in steady-state conditions 

4.1.1 Analysis of velocity fluctuations 









homogeneo 


t to a uniform static electric field E in steady 
Ltions. Diffusion and noise are related t 

lastic velocity f luctuations/v (t) of each particl 
drift value . The mathematical quantity 

ribes the common origin of diffusion and noise 
:orrelation function of velocity fluctuations, whi 
iimension, is defined as 

a t)= , 


i the brackets indicate ensemble average, and the 
i, in steady- state conditions, is independent i 
quantity carries the information on how large 
:uations are and how they decay in time. 

C (t) is related to the diffusion coefficient D tl 
:quation/4.16/: 

D = / c£fc)dt , 

o 

D can be determined by Eq.(4.2) from the evaluati< 
It is worth noting that the calculation of C(t) 
est in itself, because its analvsis. as we shall 
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Another important relation exists between the diffusion 
coefficient and the noise spectrum of velocity fluctuations, 
defined as 


U too ' L 


^ ffali) & 


U)t 

6dh 


ft-3) 


Using the Wiener-Kintchine theoremA-l?/the well known relation 


ftft) 

is found. In order to describe the different origins 

of the various terms which can contribute to the diffusion 
process of carriers in semiconductors we shall consider a 
many-valley semiconductor with two types of valleys (this is 
the case, for example, of n-Si with the external field along 
a <100 > direction). 

Let us consider an electron that, at time t, is in a 
valley of type V(t) (V(t)=l or 2) with energy 
between (r and . We may then define: 

v^=drift velocity,i,e. mean velocity of all 
electrons; 

Vy(t) =valley drift velocity,i,e. mean velocity of 
all electrons in valley V(t); 
v_y (t)=mean velocity of electrons in valley V with 
energy between £ and &A6. . 















written as the drift velocity plus a number of 
uating terms/4.12/: 

Ar/O = /fuf + jV^/ It] - lt)J + 

+ [V(t)- 'O’eWtT) - 

; nu t t y ^-5) 

^v v is the fluctuation associated with the drift 
ity of the valley in which the electron is at time 


and /to. is the fluctuation associated with the 


sing the expression in Eq.(4.1) the 
Lation function becomes: 

It) = X-I < i/\fc ((:')• otfj ( t'+1) / 


Sh c,ffc) 

■'4 J 












in the literature that the total no 













in steady-state conditions can be easily performed within 
the Monte Carlo simulation as follows. Let T be the time 
interval in which the autocorrelation function is to be 
sampled, which is usually taken as larger than the 
autocorrelation time (i.e. such that C(t)/*J0 for t^T) . 
Then T is divided into a number M of intervals of 
duration AT=T/M in order to determine C(t) at the times 

0/ AT , SAT ,MAT -T . 

During the simulation, the velocity of the sampled 
particle is recorded at the time values i-Z^T, i=0,l,2,... 
When i becomes greater than M, the products 

/ 0 r ^AT) - = fas] 

are evaluated for each i and j. Products corresponding to 
the same value of j are averaged over the simulation, thus 
obtaining 


ArCt) /vYfc+iAT) ■ C6 'at) + ^ 








The noise spectrum can be easily obtained as a Fou 
transform of C(t)/4.2/. 

In order to determine the diffusion coefficient 
means of a Monte Carlo simulation, the second central mo 


Met) 


r ~l a 

= < 2ft) - J > 


(u 


can be evaluated as a function of the simulation time 
where the average is performed over many diffe 
particles. For times larger than the initial transient, 
time dependence of M(t) becomes linear and yields 
diffusion coefficient as 


X) = 4. A- Ml(fc) 


(U 


D can also be obtained from a numerical evaluation of 
integral in Eq.(4.2) once C(t) has been obtained. 

The physical models used in the calculations for Si 
Ga As are those reported by Brunetti et al./4.21/ 
Ruch/4.22/, respectively. 


4.1.3 Results 

Results have been obtained both for covalent (si! 
case) and polar (GaAs case) materials. In order to sim| 
the complexity of the interpretation of these results 1< 








electrons in Si"have been obtained with an 


applied electric field E=10kV/cm and a crystal temperature 
of 77 K. The simplest physical situation to discuss is 
that with the valleys symmetrically oriented with respect to 
the direction of the external field (E//<?111>). 

Figs.4.1 and 4.2 report the autocorrelation function of 
velocity fluctuations and the different contributions, 
respectively, as analyzed in Sect.4.1.1, for the case E// 
^111> . Within this analysis it is seen that a negative part 









In order to connect this sequence of events to 


energy fluctuations to velocity fluctuations 6»vg. 

Fig.4.3 shows as a function of energy in Si for 

same temperature and field considered above. It can be 


here that v^ is an increasing function of £ , so tha 
€V * 

positive ££ will correspond to a positive . 


By collecting the above considerations, it can 
understood why C ^ is positive,C is negative, a 
is positive at smaller t and negative at larger t, 
a minimum which is reached at times greater than the ext 


of C 




(see Fig.4.2). 


In this case there 


the fact that the scatterirg probability is an increa 
function of energy yields a negative contribution 
longitudinal diffusivity through a negative part in 
not through a negative convective contribution , w 

exibit a regular behaviour with a small negative p 
Other off-diagonal terms are, in this case, much smaller 
Figs.4.1 and 4.4 report the autocorrelation functio 
velocity fluctuations and its different contributi 
respectively, as analyzed in Sect.4.1.1 for E// ^100)> . 
this case the different valleys have non- equiva 
orientations with respect to the field direction, 
therefore different components of the drift velocity a 
the direction of E, so that the phenomenon of longitud 
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intervalley diffusion occurs/4.21/. The intervalley 
contribution to C(t) is responsible for the long tail of the 
total autocorrelation function, absent for the case E//<flll> 
(see Fig.4.1), since the intervalley transition time is the 
largest of the characteristic times of the process under 
investigation in these conditions of field and temperature. 
If we evaluate the integral of the intervalley contribution 
C ,,we find for the intervalley diffusion coefficient the 

Qj p 

value D =24cm /sec ,which is very close to the difference D 

J IOC 

-D^=21cm^/sec /4.21/. 

The conclusion that such a difference is due to 
intervalley fluctuations is confirmed by the observation 
that the thermal and convective contributions to the total 
autocorrelation function are very similar for the two field 
orientations (see Figs<|.2 and^). 

In a transverse d irection ,for E//^100^> velocity 
fluctuations are due only to thermal fluctuations (v =v =v , 

>/ 6.V c4 

=0); their autocorrelation is always positive, since no 


enerav transfer is associated with velocitv fluctuations 
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Fig.4.5. The white-noise value of the total spectrum/ 
corresponding to a diffusion coefficient of 41 cm/sec , is 

strongly influenced by the large intervalley contribution (D w 
n 

=31 cm /sec 1. This term shows the most rapid decrease at 


the strong oscillation of C^. The convective contributi 
4 CL 

with a white-noise level corresponding to D^=10cm /sec, 
present with a monotonically decreasing behaviour, and 
always positive. The off-diagonal terms and S fc ^ have 

similar shapes of opposite signs; their cumulative 
contribution, which is relatively small, is negative at 
frequencies, corresponding to a negative contribution 
diffusivity, and becomes positive at high frequencies. 

The total noise spectrum corresponding to the sum of 
different terms seen above shows a non-Lorentian behavi 


term followed by a bump due to the thermal 
contribution and by the final dependence. Due 
to the high frequencies involved, it may be difficult to 
detect experimentally the maximum of S (U) at 800 GHz. In 


however , it must be taken into account that an 








intervalley relaxation ti 
th the above interpretat 


alculations performed 


E=200 V/cm along a^l00>and a<lll)d 


been found with experimental data of Bareikis /4.23/,as 
shown in Fig.4.6. 


Results for electrons in GaAs have been obtained with 
an applied electric field E=10kV/cm and a crystal 
temperature of 300 K. 

Fig.4.7-a reports the autocorrelation function of 
velocity fluctuations together with its three diagonal 
terms. The thermal fluctuations at the chosen field 
strength are much higher than for the case of Si, owing to 
the higher electron energy. 

Fig.4.7-b reports the off-diagonal terms. They are all 
of the same order of magnitude, but much smaller than the 
thermal contribution, so that in this part of the figure a 
different vertical scale has been used. 

The interpretation of the diagonal terms, as well as of 
the off-diagonal terms C^and , is similar to that given 
for the case of Si. Fig.4.8 shows, in fact, that v-- is a 



























large cancellation of the positive and 























Starting from an initial condition in which the parti 
positions and velocities are totally uncorrelated, 
process which occurs during the time necessary for sett 
up the correlations will be defined as the correlat 
transient. Furthermore, when a high electric field 
applied at a certain time to the electron ensemble, 
transport process itself must pass through a trans 
regime which is necessary for attaining the statio 
distribution f(^) in it-space. This process will be ca 
transport transient. 

In what follows we shall discuss how these 
different transients can be analysed separately, 
simultaneously, and how their effects influence 
transient of the diffusivity of the electron ensemble. 

The definition of the transient diffusion coeffic 
has been given by a generalization of Eq.(4.11) to arbit 
small times /4.25-4.27/: 

)*> 


where z(t) is the space position of a carrier at t 












terms of the autocorrelation function, which is easily 
interpretable from a physical point of view. By using 


Eft) 


* Zfo) > j, 


Aim clt 


(t,M) 


in Eq.(4.12)we have: 


JD(fc) • i i hi 2 W t A < fcfo) JflWt 1 )dfc' > I- 

^ t ° 

+ jdt 1 jdt" <fSirff).k/t”) > * 

oo fc 

= <kro) kft)>t Jdt' <6W $m> Mh) 

o 

If there is no correlation between the initial positions and 


velocities of the particles, Eq.(4.14) becomes: 

r k 


M)‘ JteCJz) 


fh It) 


CJz) ■- zU(D- Uli-z)>, Oi T-i L (trfi) 

where we have put in Eq.(4.15): £-=t'-t. . Eq.(4.15) reduces 

to Eq . (4.2) in steady- state conditions (t^^d)). 


Bv comDarison of these two expressions, we see that in 









inction at the t: 
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the randomized steady distribution of velocities will be set 
up, and both < v and D reach the steady-state value; the 
steady diffusivity is lower than the overshoot value because 
each particle becomes, in steady-state 
conditions alternatively slow and fast, with a reduction of 
the spreading rate at long times,as indicated by the 
negative part in the stationary autocorrelation function. 

As previously noticed, the results on transient 
diffusion can be analysed and understood in terms of the 
transient autocorrelation function C ^ ( £r) . The examples 
presented below (Figs.4.11 and 4.12) refer to the silicon 
case with E=10 kV/cm. However here the field has been 
chosen along a ^100^> direction in order to add in the 
discussion the effect of the intervalley contribution to the 
diffusivity. 

Fig.4.11 shows the transient autocorrelation function, 
as a function of fat various times t; the same initial 
conditions used in Fig...10 have been taken. Each of these 
curves is interrupted at£^t when the correlation with the 
initial conditions is reached. The area under each curve 
gives the corresponding value at time t of the transient 
diffusion coefficient. At very short times this area is 
very small both because starts from low values, and 
because it is interrupted at shortT. The first effect is a 


consequence of the transport transient , the electrons being 









still "cold", while the second effect is present because at 
small times the correlations are not yet established 
(correlation transient). As t increases, the area reaches 
larger values, and it is maximum when the positive part of C^_ 
(C) is totally present, and a negative part is not yet 
present. This leads to the overshoot of the transient 
diffusivity. When negative correlations are established, 

D(t) decreases toward the steady-state value, which is 
attained when the shape of C^(fr) reproduces the steady-state 
function and the integral of the autocorrelation function is 
extended to all significant values of . 

It is particularly interesting to reproduce the same 
analysis for the intervalley contribution alone. Fig.4.12 
shows the transient autocorrelation function for the 

intervalley velocity fluctuations (t) , defined in 

Sect.4.2.1, as a function of for various times t. The 
most striking aspect of this set of curves is the presence 
of a bump which is shifted toward a greater value of<?, as t 
increases. This phenomenon is clearly due to the velocity 
overshoot in the two types of valleys which increases the 
value of dW/t) with respect to the steady state. As t 
increases, this overshoot recedes more and more in the past, 
so that it is seen at larger correlation times until it 
vanishes, when any memory of the overshoot effect is lost. 

The overshoot of the GaAs diffusivity /4.29/can be 










ilysed in a similar way, and a negative part of the 

insient autocorrelation function prevails over the 
sitive one at the times when D is found negative/4.29/, 
s strong negative correlation is a consequence of the 
ictron transfer back and forth from central to upper 

leys, as discussed above. 

As a final result, we report in Fig.4.13 an analysis of 
transient diffusivity in a case in which the electric 
Id is applied after the onset of the zero-field 
relations. Before application of the field, the 

'fusivity slowly reaches the steady-state value in about 
psec. When the electric field is applied, a transport 
insient occurs, with heating of the carriers, on a shorter 
ie scale(in about 2 psec); during this time new 

relations are established and the diffusivity reaches the 
r lower steady value, passing through a minimum with a 
ion of negative values. This negative region is in great 
itrast with the overshoot of D seen for the case in 
.4.10. This can be explained by considering that the 
st electrons which undergo intervalley scattering after 
field application are the fastest electrons in the 
ection of the field; during the first part of the process 
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relaxation time on carrier concentration in warm electron 
conditions is observed at rather low carrier concentrations 
(10*^ -lO 1 ^ cm-3 in different semiconductors) 
/4.31-4.33/. 

The experimental data are obtained when the electron 
concentration is changed in two different ways: namely, by 
doping and compensating. In the latter case the electron 
concentration decreases when that of the impurities 
increases. Therefore the effect cannot be ascribed to 
impurity scattering, and consequently it is thought to be 
due to the interaction among carriers. For a detailed 
discussion of the experimental data see Refs.4.34,4.35. 

It is well known that e-e scattering, which is energy 
and momentum conserving, cannot have a direct influence on 
the physical mean quantities relative to the carrier 
ensemble. However the presence of carrier interaction can 
alter the probability of occurrence of other lattice 
scatterings, especially those with a well defined threshold 
of activation (e.g. emission of an optical or intervalley 
phonon)/4.36,4.37/. 

A previous comparison between experimental and 
theoretical results was performed in Ref.4.37 by means of a 






simulation of a carrier ensembl* 
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(4.17) 


















4.3.2 The Monte Carlo Procedure 



Ref. 4.42; soi 














dered. 
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(4.20) 


where g = k^-k is the wavevector difference of the collidin 
particles, k and kc are the wavevectors of the test carrie 
and the counterpart electron before their collision. 

The total scattering rate for an electron with wavevecto 
is given by 


& (£) - 


au e, I 4 T 
ziir-K 3 efo 



fOOJit c 


(4.21) 


Since the distribution in k-space of the counterpart! 
f(k c ) is not known a priori, the following rejectioi 


technique has been used. A self scattering internal to 


e-e interaction is introduced by substituting g/(g + p> ) 


Eq.(4.21) with its maximum value l/2^> 
-*> sc if * 

(k) + P (k) results to be: 


In this way 









t>e. 


sir p 3 


This new e-e total scattering probability is introduced 
among the other mechanisms. When an electron attempts an 
intercarrier collision, a counterpart carrier is selected at 
random from the particle distribution. If g is the 
wavevector difference between t!.e Monte Carlo electron and 
the counterpart, and a random number r results to be less 
than /(g**2+ ^>**2), the e-e collision is accepted. In 


this cast 


the final states of both electrons involved into 
















simulation time of the Monte Carlo electron, but thi 
difference is not relevant if the time-step duration A. 
is smaller than the time between two interparticL 
collisions• 

4.3.3 Results 

Figs.4.14 and 4.15 show typical results obtained with thi 
Monte Carlo simulation following the two approactie: 
indicated in Sect.4.3.1. 

Fig4.16 shows the dependence of ^ on the electric fiel< 
at a fixed concentration. At the highest field value: 
considered here the condition of warm-electron regime ma; 
not be fulfilled, but these fields have been considered her< 
only to check the consistence of the physical picture. 

It can be seen here that the effect of e-e interaction is 
to reduce . 

This effect can be explained through the genera! 
considerations reported in Sect.4.3.1. At low fields fe\ 
carriers have enough energy to emit optical phonons, due t( 
the low heating involved, while the majority of carriers 


interact 


ith acoustic phonons in a nearly elastic way. 
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Intercarrier collisions change the effectiveness of the 
interaction with the lattice. In particular, an essential 
contribution to the energy losses comes from two passive 
carriers (that is having energy less than A^CtJo ) if, 
during their collision, one carrier gains enough energy for 
a subsequent emisson of an optical phonon. 

As a general confirmation of this picture. Fig.4.17 shows 
the energy dissipation rate , normalised to the value 
without e-e interaction, for the optical phonon scattering 

as a function of the electric field at a carrier 

• 5 

concentration of 1.x 10 cm-3. It can be seen that at 
low electric fields the rate of dissipation is substantially 
enhanced with respect to the equilibrium case, due to the 
increase of the "active" carriers, as effect of the e-e 
scattering. At the highest fields, due to the heating, more 
carriers are active also without e-e scattering, and the 
quantities in Figs.4.16 and 4.17 become practically 
insensitive to the interaction. 

Finally in Fig.4.18 we report the normalised energy 
relaxation time as a function of carrier concentration 
compared with the experimental data of Ref.^31. 

The two theoretical techniques do not show relevant 







discrepancies in the normalised values. Even though th 
theoretical results are affected by a certain statistica 
error, there is no doubt that the effect theoreticall 
predicted is less than the one observed in the experiments 
as previously found. 

E-e interaction, even trated as two-particle collision 
cannot account for the experimental enhancement of energ 
relaxation, and more sophisticated mechanisms must bi 


invoked. 
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FIGURE CAPTIONS 

Fig.4.1: velocity autocorrelation function for electrons in 
Si, as obtained from a Monte Carlo simulation, at 77 K for 
E=10kv/cm applied along a ^100? (continuous line) and (ill) 
(dashed line) directions, respectively 

Fig.4.2: a)autocorrelation function of thermal and 

convective fluctuations for electrons in Si at 77 K and E 
=10kv/cm (see Sect.4.1.1); b) off-diagonal terms C^^ and 
Cj,£ which contribute significantly to the total velocity 
autocorrelation function for the case shown in a). 

Fig,4.3: mean velocity (continuous line and left scale) and 
distribution function (dashed line and right scale) , as 
functions of energy for electrons in Si at the indicated 
temperature and field. 

Fig,4.4: a) autocorrelation function of thermal,convective 
and intervalley fluctuations for electrons in Si at 77 K and 
E =10kV/cm (see Sect.4.1.1); b) off-diagonal terms and 

C ^ which contribute to the total autocorrelation function 
for the case shown a). 


Fig.4.5: spectral density of velocity fluctuations and its 



components for electrons in Si at the indicated temperature 
and field obtained as Fourier transforms of the 
autocorrelation function shown in Fig.4.1 and 4.4. The 
lettering of the curves is defined in the text. 

Fig.4.6: spectral density of velocity fluctuations for 
electrons in Si at 77 K and E=200 V/cm along a 0.007 and a (lll^ 
directions. Points refer to experimental data of Bareikis 
/4.23/ and lines refer to Monte Carlo calculations. 

Fig.4.7:(a) autocorrelation function of velocity 
fluctuations and its diagonal terms,and (b) off-diagonal 
terms for the case of electrons in GaAs at the indicated 
temperature and field. In (b) a different vertical scale 
has been used. 

Fig.4.8: mean velocity (continuous lines and left 
scale) and distribution function (dashed lines and right 
scale)as functions of energy for electrons in the central 
valley in GaAs at T=300K. The numbers on the curves 
indicate the field strength in kV/cm. 


Fig.4.9: spectral density of velocity fluctuations and its 
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electrons in Si at 77 K, in a physical situation in which a 
field E=10 kV/cm is switched on at time tg, =9.57 psec 
after the initial conditions (t=0) of uncorrelated 
particles. 

Fig .4.14:Carrier heating during the transient transport with 
field E=70 V/cm and T=100 K. The continuous curves have 
been obtained by a best fit of the simulation data with 
Eq . (4.17) . Numbers indicate carrier concentrations in units 
of 10**15 cm-3. Points indicate simulation data for 
n=5xl0**15 cm-3. 

Fig .4.15:Autocorrelation function of energy fluctuations in 
steady-state conditions without (closed circles) and with 
(open circles) e-e interaction for a concentration 
n=5xl0**14 cm-3. T=100 K; E=70 V/cm. 

Fig.4.16:Relative energy relaxation time as a function of 
field strength at a carrier concentration n=lxl0**15 
cm-3. Triangles and circles refer to data obtained by 
fitting the transient heating and with the energy 
autocorrelation function in steady-state, 
respectively. T=100 K. 


Fig .4.17:Relative energy loss through opt 















field strenght. T=100 K. 


Fig.4.18:Relative energy relaxation time as a function of 
carrier concentration . Triangles and closed circles refer 
to data at a field strenght E=70 V/cm obtained by fitting 
the transient heating and with the energy autocorrelation 
function, respectively. The cross indicates values obtained 
at E=20 V/cm with the two techniques, indistinguishable in 
this scale. Open circles refer to the experimental data in 


Ref .4.31. 
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CONCLUSIONS 


The present report describes the research perfor 
he Department of Physics of the University of Mode 
he ERO contract number DAJA45/83/C/0039 "Monte 
nalysis of quantum transport and fluctuatio 

emiconductors". 

As regards the analysis of velocity and 

luctuations of charge carriers in semiconductors, t 
esults have been obtained and published (they are des 
n Chapter 4). In particular, both transien 
teady-state regimes have been studied and the eff 
oise of the various sources of fluctuations hav 
nalised. 

The effect of carrier-carrier interaction on 
luctuations and therefore onenerqy relaxation time ha 


that 

observed 

As 

regard 

study 

quantui 


of the Monte Carl 















thod to the solution of the Liouville equation for 
urn transport of electrons in a phonon field has been 


tigated. The set up of the numerical procedure for this 


of practical importance is under development and 








